Localized charge bifurcation in the coupled quantum dots 
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We theoretically analyzed localized charge relaxation in a double quantum dot (QD) system 
coupled with continuous spectrum states in the presence of Coulomb interaction between electrons 
within a dot. We have found that for a wide range of the system parameters charge relaxation 
occurs through two stable regimes with significantly different relaxation rates. A certain instant of 
time exists in the system at which rapid switching between stable regimes takes place. We consider 
this phenomenon to be applicable for creation of active elements in nano-electronics based on the 
fast transition effect between two stable states. 
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I. INTRODUCTION 

Nano-scale electronics is currently a very active area of 
research. One of the main goals in this field is to design 
and characterize low dimensional structures that could 
be active elements in electronic circuitry Sin- 
gle semiconductor QDs which are referred as "artificial" 
atoms 3,0 and coupled QDs - "artificial" molecules 
, @ are usually suggested as perspective structures that 
may serve for creation of extremely small devices because 
of the possibility for the electrons spatial confinement 
on a scale less than 10 nm due to the growth processes 
01 • Double QDs systems behaviour intrinsically differs 
from a single QDs because of the variable interdot tun- 
nehng coupling which is the physical reason for 

non-linearity formation and consequently for existence of 
such phenomena as bifurcations Q and bistability [lo| . 
That's why double QDs can be applied for logic gates 
fabrication based on the effect of ultrafast switching be- 
tween intrinsic stable states. During the last decade ver- 
tically aligned QDs have been fabricated and widely stud- 
ied with the great success (for ex amp le indium arsenide 
QDs in gallium arsenide) [Tlj . [T2j . [l3|. The applied gate 
voltage dictates the electron occupancy of each QD via 
a nearby electron reservoir and tunes the relative energy 
separation between the electronic states of the two QDs 
[ll| . It was demonstrated experimentally that the local- 
ized states with different charge and spin configuration 
can be tuned into the resonance with external optical 
field due to the presence of Coulomb interaction within 
QDs These effects can lead to inverse occupation 

of different localized states, so, fully controllable solid- 
state single-emitter laser can be produced on the system 
of coupled QDs |13j. Lateral QDs seems to be a bet- 
ters candidate for scaling up the electronic coupling from 



two or several QDs by applying individual lateral gates. 
That's why they are intensively studied in the last sev- 
eral years ho\h experimentally and theoretically [l^.flH. 
Nano-devices that exhibit fast switching are supposed to 
be a basis of future oscillators, amplifiers and other im- 
portant circuit elements. The technological problem of 
QDs integration in a little quantum circuits deals with 
the careful analysis of non-equilibrium charge distribu- 
tion, relaxation processes and non-stationary effects in- 
fluence on the electron transport through the system of 
QDs [III, [131, El, El- Electron transport in such sys- 
tems is strongly governed by the presence of Coulomb 
interaction between electrons within a dot and of course 
by the ratio between the QDs coupling and coupling to 
the leads [20|. So the problem of charge relaxation due 
to the tunneling between QDs coupled with continuous 
spectrum states in the presence of Coulomb interaction 
is really vital. 

Intrinsic bistabilities in different tunneling structures 
were widely studied experimentally and theoretically. 
Obtained results provide evidence for various molecular 
[l|,[il and QDs [H , [Hi , [H , [li] switching effects appar- 
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ent in I-V characteristics due to the bistabilities in the 
tunneling current flowing throirgh the system. It was ex- 
perimentally demonstrated [ij , [^ that switching strongly 
depends on the choice of contacts, substrates and can be 
observed even for simple molecules. Coulomb interac- 
tion in such systems can be a reason for transitions be- 
tween stable states. The role of Coulomb interaction in 
double QDs bistability formation was experimentally in- 
vestigated in [l3|. Authors revealed that double-barrier 
resonant-tunneling structures have intrinsic bistable be- 
haviour in I-V characteristics due to nonlinearities intro- 
duced by the Coulomb interaction and demonstrate two 
branches with high and low current for the same voltage. 
Theoretical investigations of bistable behaviour in tun- 
neling structures usually deal with slave-boson technique 
[1^], .25], 26], drift-diffusion approach [233, [27] or Hub- 
bard approximation with negative values of Coulomb 
interaction U 12911. 
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The further progress in electronics will depend upon 
understanding intrinsic mechanisms for molecules and 
coupled QDs reversible switching from low to high cur- 
rent states but this question is still not well under- 
stood. Moreover this effect can't be accurately con- 
trolled. That's why we consider bifurcations in the sys- 
tem of QDs to be much more suitable mechanism for 
fast switching circuits creation. Bifurcation means that 
a system has several stable states or evolution regimes 
separated in time. Conditions for ultrafast switching be- 
tween these states can be controled by means of changing 
energy levels positions in QDs, the value of Coulomb in- 
teraction and strength of QDs coupling. 

In this paper we consider charge relaxation within cou- 
pled QDs due to the tunneling to continuous spectrum 
states in the presence of Coulomb interaction between 
electrons within a QD by means of Keldysh diagram 
technique |30]. Tunneling to the continuum is possible 
only from one of the QDs. We have found that for a 
wide range of the system parameters charge relaxation 
occurs through two stable regimes with different relax- 
ation rates. At a certain instant of time system switches 
rapidly between the regimes. 



II. THEORETICAL MODEL 

The model under investigation deals with the system of 
two coupled QDs with energy levels Si and £2 correspond- 
ingly (FiglD). QD with energy level £2 is also connected 
with the continuous spectrum states. Hamiltonian of the 
system can be written as: 

H = eiclci + e2c\c2 + ^ SfeC^Cfe + 

fe 

+T(c+C2+C+Ci)+^rfc(c+C2+C+Cfc) (1) 

fe 

where T and arc tunneling transfer amplitudes be- 
tween the QDs and between the second QD and continu- 
ous spectrum states correspondingly. By considering the 
constant density of state in the continuous spectrum 
the tunneling relaxation rate 7 is defined as 7 = 7rj/°T|. 
Ci"/ci(cJ/c2) and /ck- electrons creation/annihilation 
operators in the first (second) QD localized state and in 
the continuous spectrum states {k) correspondingly. We 
also take into account on-site Coulomb repulsion in the 
quantum dot with energy level £1 (first QD). Interaction 
Hamiltonian has the form: 



Hint = C/?11ct"1- 



(2) 



Different relaxation regimes are determined by the re- 
lations between the model parameters: T, 7, U and e^. 
QDs electronic states coupling T is determined by the dis- 
tance between the QDs, the barrier height and localized 
states energy levels positions. Energy levels positions are 
strongly connected with the QDs geometry: the width 
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FIG. 1: Schematic diagram of energy levels in the system of 
two coupled QDs. Second QD is also connected with contin- 
uous spectrum states. 



and the depth of potential well associated with each dot. 
The Coulomb interaction between localized electrons is 
governed by the localization radius of electronic states 
within the QDs. The value of tunneling rate 7 depends 
on the width and height of the barrier which separates 
second QD from the lead. So all these parameters can be 
varied in the real experimental situation. Our model pa- 
rameters correspond to the experimental situation when 
one of the vertically stacked interacting QDs is deep and 
narrow (deep energy levels and Coulomb interaction must 
be taken into account) and another one is wide and shal- 
low (shallow energy levels and Coulomb interaction is 
small and can be neglected) [3l| . 

As we are interested in the specific features of the non- 
stationary time evolution of the initially localized charge 
within the coupled QDs, we'll consider the situation when 
condition {si — eF)H) » 1 is fulfilled. It means that 
initial energy levels are situated well above the Fermi 
level and stationary occupation numbers in the second 
QD in the absence of coupling between the QDs is of 
order ')/{£2~'£f) « 1- So the Kondo effect can't appear 
in the proposed model and we can also omit the terms 
corresponding to the stationary solution. 

First of all we shall find localized charge relaxation laws 
in the coupled QDs without Coulomb interaction between 
localized electrons in the first QD. Let us assume that at 
the initial moment all charge density in the system is 
localized in the first QD and has the value ni(0). In the 
absence of tunneling between the QDs Green functions 
Gfi{t — t ) and G^^it — t ) can be found from expressions: 



Gfi(t-t') 



-ie(t-i')e-*"i(*-* ) 

-^e(^-^')e-*^^(*-*')-^(*-*') (3) 



where 7 ~ my^T^ is tunneling relaxation rate from the 
second QD to the continuous spectrum states. 

Retarded electron Green's function Gfi determine 
spectrum re-normalization due to tunneling processes be- 
tween QDs and can be found exactly from the integral 
equation: 



•-^11 — "-^11 +^11-' '-^22'-^ll 



(4) 
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Acting by inverse operators G°f ^ and G22 ^ integral 
equation @ can be also presented in the equivalent dif- 
ferential form (except the point t — t ): 



- 62 + ei) - T^)G^,{t, t) = (5) 

Finally, retarded Green function G{\ can be written in 
the following form: 



Gfi(t,o = *e(t-0( 



El - £2 + ij 
El — E2 
E2 — £2 + il 
El — E2 



) (6) 



where eigenfrequencies £'1,2 are determined by the 
equations: 



{E - ei){E ~ £2 + ii) - T"^ 



E 



1,2 



:(ei +' 



Z7) ± -V(ei-£2+n)'+4r2 (7) 



Let us now analyze time evolution of the electron den- 
sity in the considered system which is governed by the 
Keldysh Green function Gfiit.t) [s^: 

G<i{t,t)=ini{t) (8) 
Equation for Green function Gfi has the form: 

G<i{t,t) = Gl< + G\<T'G^2G^i + 

-t^ii 1 tJ22'-^ll + "-^ll ^ '-^22'-^ll 

(9) 

and after acting by G^f ~^ can be re- written as: 



\f^ll — 1 tJ22J'-^ll — ^ '-^22'-^ll 



(10) 



Green function Gfi{t,t) is determined by the sum of 
homogeneous and inhomogeneous solutions. Inhomoge- 
neous solution of the equation can be written in the fol- 
lowing way: 



G<(i,t')=T2 f dti f dt2G^i{t-ti) 
Jo Jo 



■^'^2(^1 ~ h)Gf^i{t2 — t ) 



(11) 



If G^2(0,0) = 0, Green function Gfi{t,t) is defined 
by the solution of homogeneous equation. Homogeneous 
solution of the differential equation has the form: 



G<i{t,t)^fi{t )e-*^^* + /2(Oe 



-iE2t 



(12) 



It is also necessary to satisfy the symmetry relations 
for function G^{t,t ): 



{G<i{t,t)r = -G<i{t,t) 



(13) 



We can determine all the coefficients because the solu- 
tion has to satisfy homogeneous integro-differential equa- 
tion: 



Gi<i(i',t) =^Ae^*^l*+*^i* -K«Be-''^i*+''^^** + 

+iB*e~'^^'+'^'^'' + iCe-'^^'+'^'i'' (14) 

We also have to fulfill initial condition: 



G<(0,0) = m° 



(15) 



As far as solution has to satisfy homogeneous integro- 
differential equation, after some calculations one can find 
the following proportionality between fi{t ) and f2{t ): 



flit) 
hit') 



£2 -El 



£2 — E2 



Z7 



(16) 



Finally time dependence of filling numbers in the first 
QD ni{t) can be written as: 



ni(t) 



i"i ■ (Ae-^(^i-^i'* +2i?e(Be-^(^i-^^*)*) 



+ C'e-^(-^^-^J)*) (17) 
where coefficients A, B and G are determined as: 



A 
B 



\E2 



r;G = 



\E^ 

\E2-Ei\^' \E2 

{E2-£i)iE*i^£l) 

\E2-Ei\-^ 



Ei\ 



(18) 



Time evolution of electron density in the second QD is 
determined by the Green function G^jC^i^ ) with initial 
condition G^2(0i 0) = 0. Green function G22it, t ) can be 
found from equation similar to the equation (|10p with the 
following indexes changing (1 o 2). Due to the initial 
conditions 712(0) = 0, ni{0) — Hq — 1, filling numbers 
evolution in the second QD n2 (t) is defined by the inho- 
mogeneous part of the solution. So time dependence of 
the electron filling numbers in the second QD n2(i) can 
be written as: 

n2{t) = {De-''^^'-^'^^' + 2Re{Ee-'^^'~^^^') + 

+ ^e-*(^^~^^')*) (19) 

where coefficients D, E and F are determined by ex- 
pressions: 
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D = F = -E ^ 



(20) 



It is clearly evident that three typical time scales exist 
in the considered system in the absence of Coulomb inter- 
action between localized electrons, which are described 
by the expressions p7)) .(??). Two of them we shall iden- 
tify as the first \Ei — and second \E2 — £'2! mode 
correspondingly. One more time scale is defined by the 
expression \Ei — £^2! ^-i^d results in formation of charge 
density oscillations in both QDs, when the following ra- 
tio between T and 7 is valid: T/7 > l/\/2- Several time 
rates in localized charge relaxation in a QD coupled with 
the thermostat were also found and carefully analyzed in 

It is necessary to mention that the suggested model 
can be generalized for the situation when both QDs are 
connected with the leads. All the expressions up to the 
equation (j20[) continue being valid if one substitutes ei 
by the expression ei — iji where 71 is a tunneling rate 
from the first QD to the contact lead. In this paper 
we are interested in the specific features of charge relax- 
ation processes first of all due to the charge redistribu- 
tion between the coupled QDs. So we consider strongly 
asymmetric case when 71 << 7 and 71 « T. In such 
approximation obtained results can be applied to the sys- 
tem of coupled QDs connected with the both leads. The 
presented assumptions correspond to the well known ex- 
perimentally vertically aligned geometry of the coupled 



QDs [lll,[ia,[l3|. 

Now we shall take into account on-site Coulomb repul- 
sion within the first QD. We shall confine ourself by an- 
alyzing only paramagnetic case when ni^ = ni-„ — ni. 

Coulomb interaction within the first QD is considered 
by means of self-consistent mean field approximation (33l | . 
It means that in the final expressions for the filling num- 
bers time evolution P7|) . (|19p one have to substitute en- 
ergy level value ei by the value ei, which is determined 
as: 



ei = ei + U + U ■ {< ni{t) > -1) 



(21) 



We consider charge relaxation from initially filled elec- 
tronic state in the first QD, and it is reasonable to deter- 
mine initial energy level position in the first QD as ei+U^ 
where ei is the energy of the empty electronic state. So 
the initial detuning is Ae = ei + [/ — £2 because at t = 0: 
ei = ei + U. 

Consequently one should solve self-consistent system of 
equations d?]), ([T7]), (dHJ and (1^ to obtain the new en- 
ergy level position Ii and ni {t). First of all it is necessary 
to substitute expressions for Ei and E2 from eq. ([7]) to the 
ea.([T8| and to determine coefficients A, B and C. Than 
one have to substitute A, B and C from ea. p^ to the 
eg. (1171) . Finally two equations are obtained where equa- 
tion for ni{t) depends on the new energy level position 



and equation for £1 depends on the filling numbers time 
evolution ni{t). These two equations can be solved self- 
consistently. The result of self-consistent solution gives 
us ni{t) and new energy level position as a functions of 
time. After this procedure coefficients Z?, E, and F can 
be found from eq. ([20]) . And finally substituting eq. ipO)) 
to eq. p^ one can obtain n2{t). 



III. CALCULATION RESULTS 

We shall start our discussion from the resonant case 
when energy levels in the both QDs are close to each 
other £1 + U ~ £2. Figl2^,b demonstrates filling num- 
bers (localized charge) time evolution in the first and 
second QDs {ni{t) and n2{t)) for the different values of 
Coulomb interaction. It is evident that Coulomb inter- 
action results in the increasing of the relaxation rate 
(grey line on Figl^h'jb) in comparison with the situa- 
tion when Coulomb interaction is absent (black line on 
Figl2^,b). A particular value of Coulomb interaction ex- 
ists in the system {U = 4 for a given set of parameters). 
When Coulomb interaction is lower than [7/7 — 4, lo- 
calized charge relaxation occurs monotonically (Figl2^,b 
grey line). Otherwise one can clearly see that relax- 
ation process becomes non-monotonic and reveals several 
typical time intervals with different values of relaxation 
rates (black dashed line on Figl2^,b). Obtained results 
strongly correspond to the localized charge relaxation pe- 
culiarities in the system of coupled QDs when Coulomb 
interaction is taken into account only within the second 
QD [20]. 

Now let us analyze non-resonant case when difference 
between the energy levels is about the values of parame- 
ters T and 7. We'll consider different signs of the initial 
detuning between energy levels. If the detuning has pos- 
itive value {ei + U > £2) for the small values of Coulomb 
interaction ([/ — 0.5, U — 1.0) filling numbers relax- 
ation rate in the first QD increases in comparison with 
the case when Coulomb interaction is absent (Fig|3^). 
With the increasing of Coulomb interaction value for 
£2 > £1 the condition \U{ni{t) — 1)| = A£ is fulfilled. So, 
energy levels detuning turns to zero at the particular time 
moment. Consequently, resonant tunneling takes place 
and charge relaxation rate reaches it's maximum value 
(L//7 — 1.5 on the Fig|3^). With the further increas- 
ing of the Coulomb interaction detuning quickly turns 
to zero and changes the sign. It results in the decreas- 
ing of relaxation rate (C//7 = 2.5 on the Fig. 3a). In the 
opposite case of negative initial energy levels detuning 
{ei + U < £2) Coulomb interaction results in the increas- 
ing of the detuning value and decreasing of the filling 
numbers relaxation rate in the first QD (FiglS};). Lo- 
calized charge relaxation in this case reveals two time 
intervals with different typical relaxation rate's scales. 
Relaxation rate in the first time interval exceeds relax- 
ation rate in the second one. 

In the case of strong Coulomb interaction (C//7 = 4 
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FIG. 2: Filling numbers evolution in the first a),c). ni{t) and second b),d). n2{t) QDs. Panels a), and b). correspond to the 
resonant tunneling (ei — £2)/^ ~ 0.0. U — 0-black line, U — 2-grey line, U — 4-dashed line. Panels c). and d). correspond 
to the tunneling between QDs in the presence of detuning (ei = e2)/7 = 1.0. U = 0-black line, U = 2-grey line. Parameters 
r/7 — 0.6 and 7 = 1.0 are the same for all the figures. 
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FIG. 3: Filling numbers evolution in the first a),c). ni{t) and second b),d). n2{t) QDs. Panels a), and b). correspond to the 
positive detuning (ei — e2)/7 ~ 1.0. Panels c). and d). correspond to the negative detuning (ei — 62)/"/ = —1.0. Parameters 
r/7 = 0.6 and 7 = 1.0 are the same for all the figures. U = 0-black line, U — 0.5-grey line, U — 1.5-black dashed line, 
U = 2.5-grey dashed line. 



on the Figl2j;,d) one can distinguish three time intervals 
with different typical relaxation rate's scales in the elec- 
tron filling number time evolution. These coincides with 
the results obtained for the system of coupled QDs when 
Coulomb interaction between electrons is taken into ac- 



count within the second QD (20j . 

Let us now focus on the most significant result ob- 
tained for the system under investigation. It is clearly 
evident that when energy levels detuning strongly ex- 
ceeds Coulomb interaction, localized charge relaxation 
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occurs monotonically with a single typical value of re- 
laxation rate (black dashed line on Fig|3^). We revealed 
that when the Coulomb interaction value within the first 
QD exceeds the values of tunneling transfer rates (T and 
7) and becomes equal to the value of strong positive de- 
tuning between energy levels (Ae >> T, 7) in the QDs 
(Ae/7 ~ U/^) or exceeds it, charge relaxation occurs 
through two stable regimes which are characterized by 
significantly different relaxation rate's values (black and 
grey lines on FigU]). 

When above-stated conditions are valid at a certain 
instant of time charge relaxation rate changes discontin- 
uously between two stable values. This phenomenon in 
the localized charge evolution can be called bifurcation. 

With the increasing of Coulomb interaction bifurcation 
takes place for the smaller time values (black and grey 
lines on Fig|4^). So one can tune the bifurcation mo- 
ment by changing the detuning value and the strength of 
Coulomb interaction. 

We have not revealed bifurcations in the case when 
Coulomb interaction was taken into account within the 
second QD [23|. The following physical reason as an ex- 
planation of this fact can be considered: Localized charge 
relaxation in the first QD occurs only due to the tun- 
neling processes to the second QD. Relaxation from the 
second dot is possible due to the coupling between the 
QDs and also to the transitions to continuous spectrum 
states. So when Coulomb interaction is taken into ac- 
count within the first QD the charge relaxation demon- 
strates much more rough behaviour due to the only one 
relaxation channel from this QD. It results in the bifur- 
cations formation. 

For the detailed analysis of charge relaxation processes 
we shall carefully examine power law exponents evolu- 
tion, which determine charge relaxation rates changing 
in each mode of the QDs (Fig|4l;) . Moreover we shall an- 
alyze time evolution of preexponenial factors which reveal 
charge distribution among the modes (FigHJi). 

Power law exponents evolution is the same for the both 
QDs (Figlll;). When the parameters values correspond 
to the bifurcation regime in the charge evolution, relax- 
ation rates of the first and second modes demonstrate fast 



transitions between the two stable values. First mode re- 
laxation rate decreases and second mode relaxation rate 
increases. After bifurcation both modes reveal identical 
relaxation rates values (Fig|31;). 

Let us now analyze preexponential factors time evo- 
lution (mode's amplitudes) in the presence of Coulomb 
interaction. In the second QD time evolution of preexpo- 
nential factors is determined by the same law (expression 
Uni (FigHJl). Time evolution of the preexponential fac- 
tors in the first QD significantly differs (expression [T51) . 
When the Coulomb interaction value becomes equal to 
the detuning value, modes amplitudes demonstrate fast 
switching between the two stable values. First mode am- 
plitude in the first QD rapidly decreases. Second mode 
amplitude in the first QD and both mode's amplitudes 
in the second QD increase. After bifurcation all mode's 
amplitudes admit equal values. It means that Coulomb 
interaction leads to equal charge re-distribution among 
the modes in the system of coupled QDs. 



IV. CONCLUSION 

We have analyzed time evolution of localized charge in 
the system of coupled QDs in the presence of Coulomb in- 
teraction between electrons within a QD. We have found 
that Coulomb interaction strongly modifies the relax- 
ation rates and the character of localized charge time 
evolution. It was shown that several time ranges with 
considerably different relaxation rates arise in the system 
of the two coupled QDs. We revealed that at a certain in- 
stant of time the system switches rapidly between stable 
relaxation regimes in a wide range of system parameters. 
This time moment can be tuned by changing the detun- 
ing and Coulomb interaction values. We consider this 
phenomenon to be applicable for active nano-electronics 
elements creation based on the effect of fast switching 
between the two stable states. 

This work was partly supported by RFBR and Leading 
Scientific School grants. 
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